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Abstract 

The effects of rare regions on the critical properties of quantum antiferromagnets 
with hidden degrees of freedom within the renormahzation group is discussed. It 
is shown that for "constrained" systems the stabihty range on the phase diagram 
remains the same as in the mean-field theory while for "unconstrained" systems the 
stability range is effectively decreased. 
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1 Introduction 



The influence of quenched disorder on the critical properties of itinerant quan- 
tum magnets is a very interesting problem in phase transition theory. Particu- 
lar interest was given to locally ordered spatial regions ( "rare regions" ) which 
are formed in the presence of quenched disorder even when the bulk system is 
still in the paramagnetic phase [1]. The conventional theory [2] ignores these 
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rare regions. For classical magnets, it has recently been shown that the rare 
regions induce a new term in the action which breaks the replica symmetry 
and that the conventional theory is unstable with respect to this perturbation 



The problem of rare regions for the case of quantum phase transitions has been 
considered in Refs. [4,5]. These authors showed that the rare regions destroy 
the fixed point found in the conventional theory of itinerant antiferromagnets 
[6] and, in contrast, critical behavior of itinerant ferromagnets [7] is unaffected 
by the rare regions due to an effective long-range interaction between the order 
parameter fluctuations. 

The conventional theory of a random- Tc quantum antiferromagnet with the 
hidden degrees of freedom on which various constraints are imposed has been 
considered in Ref. [8] and it has been shown that for "constrained" systems 
the stability range on the phase diagram remains the same as in the mean- 
field theory while for "unconstrained" systems the stability range is effectively 
decreased. In this paper wc study the effects of rare regions on the quantum 
phase transition in antiferromagnets with hidden degrees of freedom. 



2 Renormalization group equations 

Let us consider a disordered itinerant quantum antiferromagnet in which the 
n-component vector order parameter S(r, r) is coupled with the scalar non- 
fluctuating parameter y(r, r). We use an action 
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Here the function t{r, r) is the Fourier-transform of the two-point interaction 
of a quantum antiferromagnet [9,6] 



t{ci,UJn) ^t + q^ +UJ, 



(2) 



where t denotes the distance from the quantum critical point, q is the wave 
vector and cun is the Matsubara frequency. 

As the given model is phenomenological, we suppose that the coupling con- 
stant of the order parameter with the nonfluctuating parameter is purely imag- 
inary as it occurs in the Hubbard model [10-12]. It is easy to consider the case 
with a real coupling constant. The distinction between these two cases consists 
only in the definition of the part of an appropriate phase space which corre- 
sponds to a real physical model. In contrast to the classical phase transition 
the order parameter depends on both the d- dimensional vector of space r and 
imaginary time r [9]. 

In order to write an action of the conventional theory the replica trick should 
be used and then an integration over random variables should be performed. 
For simplicity we consider the case when the coefficient of the quadratic term 
in the order parameter in the action, St{r), is the fluctuating Gaussian vari- 
able with zero mean and variance A. Finally we obtain an action which has 
homogeneous saddle points only [8]. 

To incorporate rare regions into the theory we need a different approach. In 
analogy to Refs. [3-5], we consider inhomogeneous saddle-point solutions for a 
fixed realization of disorder. The partition function can be written as the sum 
of all contributions obtained from the vicinity of each saddle-point [5]. Then 
we can average over disorder by means of the replica trick. We finally obtain 
the following effective action 
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Here a and (3 are replica indices. The x-tcrm is generated by taking into 
account the inhomogeneous saddle points. The conventional theory misses 
this term. 

Let us assume the temperature T to be equal zero and use double e expansion 
[6] according to which the space dimensionality is equal 4 — e and the dimen- 
sionality of imaginary time is e,-. Of course, for the real physical case we have 
e = Er = 1. 



Defining x — xT , and putting T = 0, we obtain the following renormaliza- 
tion group flow equations in the one-loop approximation 



du 

Id 
dA 

~df 

dx 

U 

dz 

Id 
dw 

"dl 



-■ (e - 2er)u - 4{n + 8)u^ + 6uA, (4) 

: eA - 8(n + 2)?/A + 4A^ + 8nAx, (5) 

: (e - 2er)x - 8{n + 2)ux + Anx^ 6Ax, (6) 

: (e - 2er)z - 8{n + 2)uz - 2nz^ + 2zA, (7) 

: (e - 2er)w - 8(n + 2)uw - 2nw^ - Anzw + 2wA, (8) 



where I = Inb with b the scale parameter, and we have scaled K^u 
UjK^A — i> A, XjK^z —>■ ZjK^w —>■ w with = I/Stt^. We also 

denote here 



« = A.-|. (9) 

In equations (4)-(8) we separate the coefficient of the nonfiuctuating parameter 
?/(q = 0) from y{q ^ 0) because of its possible role in constraining systems 
[10,13]. 
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3 Fixed points 



Before discussing the renormalization group analysis we consider the mean- 
field theory result. After integrating over the nonfiuctuating parameter in the 

equation for the partition function wc can obtain a new effective action in 
terms of the order parameter S(r). In the mean-field approximation for this 
effective action the boundaries of a stability range can be easily found. It 
is convenient to introduce new notations for the coupling constants of the 
effective action as 

Ai°) = A«=/.V2/?. (12) 

For the unconstrained system A^^-* = A^^-* = Ac (or z = 0) while the boundaries 
of the stability range correspond to the equations of lines 

A, - A(°) = 0, Af = 0. (13) 

For the constrained system A^°) = (or w = 0) while the boundaries of the 
stability range can be written as 

A, = 0, AW = 0. (14) 

The stability ranges are represented on planes (Ac — A^) (Fig. 1) and (A^^^ — A^) 
(Fig. 2). 

The characteristic feature of fiow equations is the closed system of three equa- 
tions (4)- (6). Within the notations this set of equations coincides with the 
appropriate set of equations of Ref. [4] where however the nonfiuctuating de- 
grees of freedom were not taken into account. It is easy to find all eight fixed 
points (4)-(6). Four of the fixed points have a zero fixed-point value of x*, 
X* — 0. The other four fixed points have x* ^ 0. 

Using these fixed points we can find other fixed points from the fiow equations 
(7), (8). There are the sixteen fixed points with x* = (Table 1) studied before 
in Ref. [8]. The other sixteen fixed points have x* ^ (Table 2). As can be seen 
from Table 1 we can separate all fixed points on groups consisting of two fixed 
points, for example group of Gaussian points (G) and group of renormalized 
Gaussian points (RG). Each group has its proper set of critical exponents. 

Such classification of fixed points was used for description of phase transition 
in classical systems [10,14,15] where for renormalized fixed points the criti- 
cal exponents can be found according to Fisher [16] via critical exponents of 
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Table 1 

Fixed points with x* = 



u* A* X* z* w* 



e-26r 

2n 



Gi 

G2 



RGi 

RG2 ^ 

Hi 

H2 ^ 

^ ^ 2n(n+8) ^ 

RHn g-2er n n (e-26^)(4-n) {e-2er){4-n) 

4(n+8) ^ ^ 2n(n+8) 2n(n+8) 

Ui _| ^ 

U2 -| 

RUi -I ^ 

RU2 -f ^ 

p e+4eT (4^n)£+4(n+2)eT- {n-4)e-12ner 

■^1 16(n-l) 8(n-l) " ^ 8n(n-l) 

Rn ^+4£r (4-n)£+4(n+2)£^ H 

^2 16(n-l) 8(n-l) U U U 

RR e+4eT (4-w)e+4(n+2)e.r (n-4)e-12ner 

16(n-l) 8(n-l) ^ 8n(n-l) ^ 

RR e+46T (4— n)e+4(n+2)£T- p, (n— 4)£— 12neT- (n-4)e-12neT- 

-'-^-'-^2 I6(n-1) ' 8(n-l) ^ 8n(n-l) 8w(w-l) 

appropriate non-renormalized values 

'-'^renorm — i ^renorm ^ • 

1 — a 1 — a 



The eigenvalues Aj (i = 1, . . . , 5) of flow equations linearized about the fixed 
point (4)- (8) are indicated in Table 3 and Table 4. 

It should be noted that eigenvalues Ai and A2 for random (R) and renormalized 
random (RR) fixed points for n > 1 are complex. 

Due to the relation 

ii = A,-A«, (16) 
the fixed points with A* = and x* — ^ ahgn on parallel fines: 



(1) As = Ac for the unconstrained system (z = 0); 

(2) A5 = A^^^ for the constrained system (w = 0); 
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Table 2 

Fixed points with x* 7^ 





7;* 


A* 


r* 




Z* 


7/?* 


1 


n 

u 





e-2er 
4n 







e-2eT 
2n 


2 


n 





e-2eT- 
4n 




u 





Q 









e-2er 
4n 




e-2eT- 
2n 


u 


4 





n 
u 


e-2er 
4n 




e-2er 
2n 


e-2er 
2n 


5 


e-2tr- 





(n-4)(e-2e^) 







(e-2e^)(4-n) 


4(n+8) 


4n(n+8) 




2n(n+8) 


6 


e-2tT 


n 
u 


(n-4){e-2e^) 




n 

u 





4(n+8) 


4n(n+8) 




7 


€-2eT 


u 


(n-4)(e-2er) 




(e-2er)(4-n) 


n 

u 


4(n+8) 


4n(n+8) 




2n(n+8) 


8 







(n-4)(e-2er) 




(e-2er)(4-n) 


(e-2er)(4-n) 




4(n+8) 


4n(n+8) 




2n(n+8) 


2n(n+8) 


9 


n 
u 


-e+4e^ 


e+4er 







3e-4eT- 


8 


16n 




8n 


10 


n 

u 


-e+4er 
8 


e+4er 
16n 




u 





11 


n 

u 


-e+4er 


e+4e^ 




3e-4er 


u 


8 


16n 




8n 


12 





-e+4er 


e+4er 




3e-4er 


3e-4eT 


8 


16n 




8n 


8n 


13 


e+4er 


{n-4)e+24tr 


(n-4)(e+4eT) 







(4-n)3e-(n+8)4eT- 


8(10-n) 


4(10-n) 


8n(10-n) 




4n(10-ra) 


14 


e+4e^ 


(ri-4)e+24eT 


(n-4)(e+4er) 










8(10-n) 


4(10-n) 


8n(10-n) 




15 




(n-4)e+24er 


(n-4)(e+46^) 


(4. 


-n)3e-(n+8)4eT- 





8(10-n) 


4(10-n) 


8n(10-n) 




4n(10-n) 


16 


e+4er 


(n-4)e+24eT 


(n-4)(e+4er) 


(4- 


-n)3e-(n+8)4eT 


(4-n)3e-(n+8)4er 


8{10-n) 


4{10-n) 


8n(10-n) 




4n(10-n) 


4n(10-n) 



(3) A5 — Ac = (e — 2er) /4:K4{n + 8) for the unconstrained system (z — 0); 

(4) Xs - = (e - 2€r)/4:K4{n + 8) for the constrained system = 0). 

It is easy to see that for the quantum phase transition {tr 7^ 0) in the case 
when the inequahty e < 2er is satisfied the fixed points lying on the two last 
lines are not situated in the stability range of the effective action (Fig. 1 and 
Fig. 2). Thus the stability range for the quantum phase transition in the con- 
strained system remains such as in the mean-field theory in contrast with the 
classical phase transition where these fixed points are situated in the stabihty 
range and cause an effective reduction of this range [11,12]. However for the 
unconstrained system there is the unstable fixed point Gi on the boundary of 
the stability range. Its critical exponents coincide with critical exponents of 
the point G2- Analogous to the classical phase transition a separatrix going 
out this point divides the stability range in two regions. One is the region of 
accessibility of the fixed point G2 and the other is the region of phase space 
where flow trajectories approach to the "invariant" line G2 — Gi. 



Correspondingly, the fixed points with A* 7^ and x* 
parallel lines: 



are situated along 
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Table 3 

Eigenvalue of fixed points with x* = (Af 2 



8(n-l) 



{-3ne + 4(n - 4)er ± [(5n 



8)2e2 + 24(16 - 12n - ^2)66^ + 48(16 - 8n - 5n2)e2]V2}) 





Ai 




A2 


A3 


A4 


As 


Gi 


e - 2er 




e 


e-2er 


e - 2er 


-e + 26^- 


G2 


e - 2er 




e 


e-2er 


e-2er 


e-2ei- 


RGi 


€ - 2er 




e 


e - 2er 


-e + 2er 


-e + 2eT- 


RG2 


e - 2er 




e 


e - 2er 


-e + 2er 


e - 2er 


Hi 
Ho 
RHi 
RH2 


-e + 2er 

-e + 2e^ 
-e + 2e^ 


(4- 


-n)e+4{n+2)eT 


(4-n)(e-2e^) 


(4-n){e-2e^) 


(4-n)(e-2e^) 


(4- 


n+8 
-n)e+4(n+2)eT 


n+8 
(4-n)(e-2e^) 


n+8 
(4-n)(e-2e^) 


n+8 
(4-n)(e-2e^) 


(4- 


n+8 
-n)e+4(n+2)er 


n+8 
(4-n)(e-2er) 


n+8 
(4-n)(e-2er) 


n+8 
{4-n)(e-2er) 


(4- 


n+8 
-n)e+4(n+2)er 


n+8 
(4-n)(e-2er) 


n+8 
(4-n)(e-2er) 


n+8 
(4-n)(e-2er) 




n+8 


n+8 


n+8 


n+8 


Ui 


-|e-2er 




— e 


-\e-2er 


|e - 2er 


-|e + 2e^ 


U2 


-ie - 2er 
2 ~ 




— e 


-\e - 2er 
2 ^ 


U - 2er 

2 ^ 


he - 2er 
2 ^ 


RUi 


-ie-2er 




— e 


-\e-2er 


-ie + 2e^ 


-^e + 2er 


RU2 


-ie- 26^ 




— e 


-\e-2er 


-ie + 2e^ 




Ri 
R2 
RRi 
RR2 


Af 
Af 
Af 
Af 




Af 
Af 
Af 
Af 


(4-n)(e+4c^) 


(n-4)e— 12neT 


(n-4)e-12neT- 




4(ra-l) 
(4-n)(e+46r) 


4(n-l) 
(n— 4)e— 12neT 


4{n-l) 
(n— 4)€— 12ner 




4(n-l) 

(4-n)(e+4er) 
4(n-l) 

(4-n)(e+4e^) 


4(n-l) 

(n— 4)e— 12neT 
4(n-l) 

(n— 4)e— 12raeT 


4(n-l) 

{n—'i)e— liner 
4(n-l) 

(ra— 4)e— 12neT 




4(n-l) 


4(n-l) 


4(n-l) 
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Table 4 

Eigenvalue of fixed points with x* j^O (A^^| = |{4e^ - 3e± [256^ - 24ee^ - 2406^] V2} 
and A^la^ = -4(1^ {(16 - n)e + 4(n - 4)6^ ± [9(8 - n)^€'^ + 24(71^ - 4n - 48)ee^ + 
16(ra^ +40ra-464)e^]V^}) 





Ai 




A2 


A3 


A4 


As 


1 


e - 2er 




-e + Ae-j- 


-e + 26^- 


e-2er 


-e + 2er 


2 


e - 2er 




-e + 4ei- 


-e + 2er 


e-2ei- 


e - 2er 


3 


e - 2er 




-e + Atr 


-e + 2er 


-e + 2er 


-e + 2er 


4 


e-2er 




-e + 4e^ 


-e + 2er 


-e + 2e^ 


e-2er 


5 


-e + 2er 


(n 


-A)c+2Aer 


(n-4)(e-2er) 


(4-n)(e-2er) 


(A-n){e-2er) 




n+8 


n+8 


n+8 


n+8 


6 


-e + 2er 


(n 


-4)e+24er 
n+8 


(n-4)(e-2e^) 
n+8 


(4-n)(e-2e^) 
n+8 


(4-n)(e-2eT) 
n+8 


7 


-e + 2er 


{n 


-4)e+24er 
n+8 


(n-4){e-2e^) 
n+8 


(4-n)(e-2e^) 
n+8 


(4-n)(e-2e^) 
n+8 


8 


-e + 2er 


(n 


-4)e+24er 
n+8 


(n-4)(e-2€T) 
n+8 


(4-n)(e-2€r) 
n+8 


(4-n)(€-2er) 
n+8 


9 


4 




3,(9) 
-^2 




3€-4£^ 

4 


3e+4e,- 
4 


10 


e+4er 
4 




'^2 


3,(9) 


3e-4£r 
4 


3e-4eT- 
4 


11 


e+4eT 
4 




3,(9) 
'^2 


3,(9) 


3e+4eT 
4 


3e+4eT 
4 


12 


e+4€T 




X(9) 


3,(9) 


3e+4€r 
4 


3e-4€r 
4 


13 


a(^=^) 




3,(13) 
^^2 


(n-4)(e+4e^) 


(l-n)3c-(n+8)4eT 


(4-n)3e-(n+8)4eT 




2(10-n) 


2(10-n) 


2(10-n) 


14 


a(^^) 




X(13) 
^^2 


(n-4)(e+4er) 


(4-n)3e-(n+8)4eT 


(4-n)3c-(n+8)4cT 




2(10-n) 


2(10-n) 


2(10-n) 


15 


a(^^) 




x{13) 
^^2 


(n-4)(e+4er) 


(4-n)3e-(n+8)4eT 


(4-n)3e-(n+8)4er 




2(10-n) 


2(10-n) 


2(10-n) 


16 


aS^^) 




3,(13) 
^^2 


(n-4)(e+4e^) 


(4-n)3e-(n+8)4er 


(4-n)3e-(n+8)4eT- 




2{10-n) 


2{10-n) 


2{10-n) 




Fig. 2. Phase diagram for the constrained system (A* = 0, x* = 0) 
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Fig. 3. Phase diagram for the unconstrained system (A* 7^ 0, x* = 0) 




RRi 1^ 



Fig. 4. Phase diagram for the constrained system (A* 7^ 0, x* = 0) 

(1) A^. — Ac = (e + Aer)/lQK4{n — 1) for the unconstrained system {z = 0); 

(2) Xg — A^^^ — {^ + ^e-r)/l^Ki{n — 1) for the constrained system [w — Qi). 

For the unconstrained system the fixed points lying on these lines for both 
the classical (e^ = 0) and the quantum phase transition are situated in the 
stabihty range (Fig. 3). Thus the accessibihty range of fixed points is also 
reduced in comparison with the mean-field result. For the constrained system 
not all points are in the range of stability (Fig. 4). 

The analysis of the renormalization group equations (4), (6) and (5) shows 
that as well as in the case of systems without hidden degrees of freedom [4,5] 
the Gaussian fixed point G2 is stable for d > A and unstable to disorder for 
d < 4. The random fixed point R2 is stable for 3ne > 4(n — 4)e^, n > 4 and 
12neT- > (n — 4)e. It is easy to see that this point is always stable for 0? < 4 
and 4 < n < Tie where ric — 16 for the particular case = e. It should be 
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noted here that in the conventional theory this fixed point is stable for d < A 
and n < 4 [6,8]. For ri < 4 the stable unphysical fixed point U2 is unaccessible 
from the physical range (according to definition the magnitude A should be 
positive) and for d > 4 it is unstable. The Heisenberg fixed point H2 for the 
pure quantum system becomes unstable because the dynamic critical exponent 
z reduces the upper critical dimensionality. There also are new fixed points 
with X* 7^ 0. As can be seen from Table 4 they are all unstable except the fixed 
point 14, which is stable for n < 4 and 3(4 — n)e — 4(n + Sje,- < 0. However, 
this fixed point is also unaccessible from physical range since the value of x* 
is negative. 

Thus, we see that the quantum character of the phase transition in systems 
with hidden degrees of freedom with constraints leads to conclusion that for 
pure systems (A* — x* — 0) the phase transition is described by the Gaus- 
sian fixed point with "classical" critical exponents of the mean-field theory 
whereas for the disordered systems (A* 7^ 0, x* ^ 0) the rare regions de- 
stroy the phase transition for n < 4 similar to systems without hidden degrees 
of freedom [4,5]. However, for 4 < n < rzc the phase transition is described by 
the random fixed point. Hidden degrees of freedom result in a decrease of the 
stability range with comparison to the result of the mean-field theory and ex- 
istence of the range of the phase space in which the fiow trajectories runaway 
from fixed points not intersecting the stability range of the mean-field. 
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